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This paper is partly of a tutorial nature. The intended purpose is to exploit the essen- 
tial two-dimensional nature of wave propagation in the earth-ionosphere waveguide. It is 
shown that, without resorting to erudite arguments in the complex plane, the usual working 
formulas of VLF mode theory may be derived directly from orthogonality considerations. 
Furthermore, the physical insight gained by the present development immediately suggests 
how the formulas may be generalized to an earth-ionosphere waveguide of nonuniform width. 

1. Introduction 

In the theory of terrestrial radio propagation, the earth is represented by a spherical body 
and the at mosphere is usually idealized by concentric spherical layers. In this sense the problem 
is formulated as a three-dimensional one. However, in subsequent developments and in the 
reduction to useful formulas, approximations are usually made which display the inherent two- 
dimensional character of the problem. Typical of these approximations are the replacement of 
the spherical wave functions by Airy functions and the Legendre functions by leading terms in 
their asymptotic expansions. 

It is the purpose of the present paper to develop the final formulas directly from a two- 
dimensional formulation. This should help dispel any doubts concerning the validity of the 
approximated three-dimensional solutions. A demonstration such as this is timely since a two- 
dimensional laboratory model has been used recently to study VLF propagation in a spherical 
waveguide. 1 At the same time we shall discuss the excitation factors, height-gain functions 
[Spies and Wait, 1963], and mode conversion coefficients [Wait, 1962a] for the two-dimensional 
model. For the purposes of this paper, the ionosphere is assumed to be equivalent to an iso- 
tropic medium. 

2. Primary Fields 

In formulating the radio problem in spherical coordinates, it is customary to choose the 
source as a radial electric dipole or a radial magnetic dipole. By an appropriate superposition, 
the fields of an arbitrary source may be determined. In a two-dimensional model it is more 
convenient to work with line sources. For example, in place of a radial electric dipole on a 
spherical surface it is convenient to employ a line dipole source. To introduce the subject 
the fields, of such a line dipole located in free space, are derived from basic principles. 

The line source consists of a uniform distribution of vertical electric dipoles. The situation 
is illustrated in figure 1 where the source current density is / amps per meter, for a strip of width 
dl, along the x axis. Because the source current has only a z component, the resulting Hertz 
vector also has only a z component, II.. Clearly, the Hertz vector of the distribution is ob- 
tained by integrating the contributions from the distribution of vertical electric dipoles; thus 
[Wait, 1959] 
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i Microwave model experiments being conducted by the author's colleagues E. Bahar and S. Maley in the Department of Electrical Engineer- 
ing at the University of Colorado. 
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Figure 1. Line dipole source and coordinate system 
employed for the discussion of primary fields. 

where k = 2r/wavelength , 

r=[(x-x') 2 +y 2 + z 2 ] l/2 , 
and eo =8.854X10" 12 . The integration with respect to x f may be readily carried out, to 
give [Wait, 1959] 

n -2^^>' (2) 

where p= (y 2J r z 2 ) 1 ' 2 and K is a modified Bessel function of order zero. The resulting magnetic 
field, which has only an x component, is obtained from 
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where /x =47rX10~ 7 . In the far zone, where &p>>l, this reduces to 
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where sin a=y/p. It is of interest to observe that the radiation pattern characterized by sin a, 
is identical to the pattern of an isolated vertical dipole. However, in the case of a line dipole 
source, the radiation field is a cylindrical wave whose amplitude varies inversely as the square 
root of distance. This is to be contrasted with the radiation field of a point electric dipole 
which is a spherical wave with an inverse distance variation of the amplitude. An obvious 
point of similarity between the line dipole and the point dipole is the common phase factor 
exp (—ikp), occurring in the radiation fields. 

The preceding rather elementary derivation is based on a prior knowledge of the expression 
for the Hertz vector for a point dipole. It is instructive to obtain the same final formula for 
H x by using a transform method. 

The starting point is to recognize that H x must satisfy the wave equation 
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Solutions are of the form exp (±i\z) exp (±uy) where u=(\ 2 —k 2 ) l/2 and X is a parameter. 
A general solution is constructed by writing 
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(6) 



where F(\) is an unknown function and where the contour of integration is along the whole 
real axis of A. The radical is to be chosen in u such that Re u^>0. Thus, the minus sign in 
the exponent is to be employed when ?/>() and the positive sign is chosen when ?/<C0. Diffi- 
culties with branch points on the contour are avoided by allowing k to have a very small but 
finite negative imaginary part (i.e., — Imag. F>0). 

The function F(\) is obtained by noting that H x must become numerically equal to —1/2 
at z=0 and as y tends to zero from positive values. Thus 

H x ] st-U{z)dl, (7a) 

where h{z) is the unit impulse function at z=0. It should also be noted that 

H t ] s+|«(g)<#. (7b) 

Employing the spectral representation of the unit impulse function, i.e., 

s(2)= i rr e ~™ dx ' (8) 

it follows from (6) and (7a) that 



F(X) = - 1 ^- (9) 



Thus 
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for y^>0. This integral in (10) may be evaluated by noting that 

K (ik P )=l r e-^ z e- uv u-\l\, (11) 

" J -co 

which may be differentiated under the integral sign with respect to y. As a result, 

H __iJcM Ki(i/cp)sma (12) 

Z7T 

where sin a=y/p and K x is the modified Bessel function of order one. In the far zone, where 
Arp»l, it follows that 

^-|£f ^' sin «- (i3) 

which is identical to (4) as it should be. 

3. Fields in the Concentric Cylindrical Waveguide 

To come to grips with the general problem, a rather simple model is chosen. The wave- 
guide region consists of two concentric cylindrical surfaces of radii a and a-\-h, respectively, 
as indicated in figure 2a. In terms of a cylindrical coordinate system (r, 6, z), these bounding 
surfaces are defined by r=a and r=a+h where the tangential fields are assumed to satisfy 
impedance-type boundary conditions. The field variations in the z direction are taken to be 
zero so that the modes in the structure may be either TM (transverse magnetic) or TE (trans- 
verse electric). Only the TM modes will be considered here since their analog in VLF radio 
propagation is of greatest practical interest. The analysis for the TE modes is almost identical. 

For the TM modes in such a cylindrical region, the magnetic field has only an axial or z 
component which is denoted H (without writing the subscript z). From Maxwell's equations, 
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Figure 2a. Concentric cylindrical model and 
conventional coordinate system (r, 0). 



it is seen that the electric field has only r and components. 

The magnetic field in some aperture plane 6=d , a+A>r>a is now assumed to be known 
and designated by H(0 , r) , being independent of z. The boundary conditions for H(0, r) are 

E =-Z g H at r=a, (14a) 

E e =Z t H at r=a+h, (14b 

where Z g and Z^ are the respective surface impedances. Within the region a+h^r^a, the 
magnetic field must satisfy the wave equation 

(?i-i+?s+ t ') H <^'= - < i6) 

Solutions are any linear combinations of the form 

H™(kr) 6 ' 

where v is a complex quantity independent of the coordinates. H™ and H™ are Hankel func- 
tions, of order v, of the first and second kind, respectively. 

Solutions which are periodic require that v is an integer. However, this is an unnecessary 
restriction since individual modes need not be periodic. Instead, the value of v will be deter- 
mined from the boundary conditions. Furthermore, without loss of generality, it is specified 
that Re z>>0, and attention is restricted to solutions which behave as e~ ivd . This will correspond 
to waves which are attenuated in the positive direction. The waves propagating in the 
negative direction are identical in form. 

If attention is restricted to large radii of curvature such that both kcT>^>l and h/a<C<Cl 
are satisfied, the Hankel functions may be approximated by their Airy function representation. 
There is a vast amount of literature on this particular subject and here only the final form is 
quoted. Thus [Wait, 1962b] 

H*m^(Q H Mt-v), (16) 
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and W\ is an Airy function. Similarly, 
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Figure 2b. Concentric cylindrical model with the 
"natural" coordinates (x, y). 



in terms of the Airy function w 2 . In terms of more conventional notation 

v>i(t) = T^[Bi(f) -iAi(t)], (18) 

and 

w 2 (t) = ^[Bi(t)+iAi(t)] } (19) 

where the Airy functions Ai and Bi have been defined and tabulated by Miller [1946]. The 
functions w x and w 2 satisfy 

^)-^(0=0, (20) 

which is known as Stokes' (or Airy's) equation. 

To facilitate the subsequent discussion, certain dimensionless parameters are introduced 
which simplify the notation. These are 

x=(ka/2)X6, Xo=(ka/2)^e 0l y=[2/(ka)]^k(r-a), y =[2/(ka)]^kh, g=-i(fca/2)«Z^ , 

2«=--i(fo/2)*Z</io, and r, Q = (moAo)^120tt. 

The parameters x and ?/, while describing horizontal and vertical distances, should not be con- 
fused with Cartesian coordinates. The quantities q and q { describe conveniently the electrical 
properties of the bounding walls of the waveguide. The situation is illustrated in figure 2b. 
The (approximate) solutions of (15) are now written as linear combinations of 

v>i(t—y) 

e~ m , 
w 2 (t—y) 

The discrete values of t (i.e., the eigenvalues), denoted by t n , are determined from the boundary 
conditions. Thus, the resultant field has the form 

H=H(x, y)= S M(*», y)e- i(x - T o )t n e - m(e - d o\ (21) 

71=1,2,3 . . . 

where 

*(*», y) = Mtn-y)+A(t n )w 2 (t n -y), (22) 

where b n and A(t n ) are undetermined coefficients. It is evident from (20) that <£ satisfies 

(|r-t+v)*(t,y)=o. (23) 

The boundary conditions given explicitly by (14a) and (14b) may be written in the form 

3£*«.v>] f _ +«*(*. °)=°. < 24) 
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and 

i *(*,»)! -gLMt,y )=0. (25) 

The latter two equations determine both the eigenvalues t n and the coefficient A(t n ) . Using (22) 
and (25), it is seen that 

where the prime indicates a derivative with respect to the argument t—y of the Airy function. 
If (24) is also to be satisfied, it follows that 



A(t n )B(t n ) = l = e- i2 * n , (27) 

where 

» 2(t)—qw 2 (t) ~l 

;(*)-2«*(oJ ( } 






Within the approximations used, the modal equation (27) for this cylindrical model is identical 
to the one developed for a concentric spherical model [Wait, 1961]. 

It should be mentioned at this stage that the surface r=a J rh can be regarded as a reference 
surface where the ratio of the tangential fields is specified. In the general case, Z x or q t may 
be a function of the eigenvalues t n . However, for VLF it is a good approximation to regard 
Z t or q_i as constants. This is equivalent to stating that the surface impedance does not depend 
on the angle of incidence. In a similar manner Z g or q can be regarded as constants. 

The orthogonality properties of the modes are now studied. We consider two sets of 
values, t n and t m , which satisfy the boundary eqs (14a) and (14b). However, for any value of 
y these must also satisfy 

£- 2 <f>n-(t n -y)<3> n =0, <f> n =$(t n ,y), (29) 

and 

^ $ m -(t m -y)$ m =0, * m =*(tm,y). (30) 

After multiplying the first of these equations by $ m and the second by $ n , they are subtracted 
from one another. Both sides of the resulting equations are then integrated with respect to 
y over the range to y . This results in 

d d To rvo 

® n du * m ~~* m dy* n \ = ( tn ~ tm ) I $n$ m dy. (31) 

In view of the boundary conditions on <£> n and $ m at y=0 and y=yo, the left-hand side of the 
preceding equation is zero. Thus, the integral on the right also vanishes if t n is not equal to 
t m - Therefore, we have the important result 

J^ Ht n ,y)*(tn,y)dy=0 if m^n. (32) 

It now follows that, if both sides of (21) are multiplied by $(t m , y ) and integrated from 
to y , 

po 

H(x ,y)$(tn,y)dy 

b n =%~ - (33) 

[Ht n ,y)fdy 

Jo 



The normalizing integral 



N n = ( VO lHt n ,y)fdy, (34) 

Jo 
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is now expressed in a more convenient form. This is accomplished by noting that 

fmt n ,y)fdy=-(t n -y)mtn,y)f+{^(tn,y)} 2 . (35) 

This can be proved by differentiating both sides with respect to y and making use of (29). 
Then using the definition of $(t n , y) in terms of Airy functions, it follows that 



and 



2? 
$(*»,0) = , u , — rrv (36) 

2% 

Ht '" y(,) = ~mn-y< > )+<LrW 2 (tn-yo)' ^ 



where use is also made of the Wronskian condition 

w l (t)w 2 (t)-w[(t)w 2 (t) = -2i, (38) 

which is valid for any value of t. Finally, on making use of (35), (36), and (37) along with 
(14a) and (14b), it is found that 

at = 4(^-g 2 ) 4(t w — y --g?) , . 

n WSn)-qw 2 {t n )r[w'Sn-y,) + (liW 2 {t n -y,)f ^ 

This can be regarded as a fairly important result. 

Equation (33) for the coefficient b n can be written in the convenient form 

H(xo t y)$(t n ,y)dy 

i "A-n JO /'Af\\ 

K= ^ ~ [*(«., o)]« -' (40) 



where 



We are now in the position to obtain the fields resulting from a line dipole source at height 
z above the earth's surface. Thus, the aperture field is specified by 

H(r,0)=8(z-z)(I/2)dl, (42) 

where z=r—a and where / and dl have the same meaning as used previously. In terms of 
the natural coordinates, the assumed field in the aperture plane is given by 

H(x ,0)=k(^ i 8(y-y) I ^dl, (43) 

A / 2 \ H A 

where y=l j—j kz. Using this result and (40), it follows that 

K y. [*(*,, 0)]»W m 

The resultant field for a dipole line source, at 0o=O, may now be written compactly in the 
form 

H(x,y)=H«»(x,y)V, (45) 

where 

F= 4^ e _ iF/4 ^^ -'*>nA n G n (y)G n (y), (46) 

y<> n=i,2,3 . . . 
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and 



H ^> ^=|SP « W |? (£) V ~ (47) 

*-<*>-ffet (49) 



The quantity H (0) , defined above, is numerically equal to the broad-side field of the line 
dipole source in free space at a distance equal to ad (e.g., compare with (4) for a=ir/2). The 
quantity 6 n (y), described as a height-gain function, is the ratio of the field of mode n at height 
y to the field of the same mode at y=0. 

The coefficient A n , which is defined by (41), is called the "excitation factor/' The physical 
significance of this factor may be evident from the identity 

-^- (50) 



f 



[Gn(y)Wy 



which may be verified by comparing the right-hand sides of (33) and (40). As may be seen, 
height-gain functions G n (y) f which increase appreciably with height, correspond to small 
excitation factors. 

The expression for V given above is identical in form to the corresponding result derived 
for a spherical model under the assumption that A/&<<1 and &&>>1. In both cases, Fmay 
be described as the ratio of the actual field to a "primary field." Furthermore, the height- 
gain function 6 n (y) and the excitation factor A n also have identical definitions in the cylindrical 
and spherical models. 

4. Extension to Waveguides of Nonuniform Width 

The development of the theory in this paper lends itself quite readily to certain generaliza- 
tions. A case of some interest is when the height of the upper boundary varies in the direction 
of propagation. In the two-dimensional model considered here, the (normalized) height 
y (x') is considered to be a function of x' , the (normalized) distance. As indicated in figure 3, 
y Q (x') varies, in a smooth manner, from y (Q) above the source to y (x) above the observer. 
The (normalized) surface impedances of the walls are also regarded as a slowly varying func- 
tion of x r . 

On physical grounds it can be expected that a waveguide with slowly varying properties 
will not differ fundamentally from a waveguide of constant cross section. In other words, at a 
distance x' , the structure of the modes is characteristic of a uniform waveguide of constant 



q.(x') 
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Figure 3. Waveguide of nonuniform width. 
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width which is equal to y (x') and with constant wall impedances equal to q(x') and q t {x f ). 
To obtain the overall behavior of the field one is led to integrate over the whole range of x f 
from to x. Thus, the appropriate form of the (complex) phase term is 



exp — i t n (x')dx' exp (—ikaS), 



where x=(ka/2) H 9 } x' = (ka/2) Va 6' , and t H (z') is the slowly varying eigenvalue. From what 
we have said above, it follows that t n (x') is a solution of 

A[t n (x')]B[t n (x')]=exp {-i2nm), (51) 

where A[t n ] and B[t n ] have the same definitions as (26) and (28), where now y , q, and q { are 
functions of x' . 

For the slowly varying waveguide, it is also necessary to employ the height-gain functions 
determined by the local width. Thus, G n {y) is to be replaced by the function 

G »«>y^W)i' (52) 

which is an obvious generalization of (48). 

The resultant expression for the field is thus given by 

H(x,y)=H™(,r,y)V, (53) 

where 

f =» e ' M4 -,£... exp [-'' J> (B#) *'] ***•«>> *>*-<». ">■ <"> 

The "excitation factor" is now obtained from 

r [y (Q)yo(x)} l < 2 , ... . 

A »~~ r rv (0) rv (x) "\ 1/2 (,00a; 

2 {Jo [G, » (0 ^ )]2 ^J W&'ViWvY 

which is an obvious generalization of (50). It is then evident that 

A n =[A«»A<?] 1/2 , (55b) 

where A£ 0) is defined by (41) if y , q_, and q_ t are replaced by y (0), q(0), and ^(0), respectively. 
Similarly, the definition of A ( n x) has the form of (41) if the corresponding quantities y Q (x), 
q(x), and q^x) are employed. 

It might be mentioned that if the curvature of the walls is sufficiently small and q { is 

essentially zero, A n will become equal to unity as in the guide of constant width [Wait, 1962b]. 
Thus, in this limiting case, it may be seen from (54) that the fields of the individual modes 
vary approximately as l/-yjy (x), apart from the exponential factors. This behavior is consistent 
with conservation of power since the outward power flow is proportional to l/yo(x) or inversely 
to the width of the waveguide. 

The form of (54) is suggestive of a WKB approximation for propagation in a horizontally 
stratified medium with slowly varying properties. Implicit in the development is the assump- 
tion that conversion of modes from one order to another may be ignored. Thus, the properties 
of the waveguides must vary only slightly in a distance equal to one wavelength. 

5. Waveguide With Localized Obstruction 

If the width of the waveguide changes suddenly, it is apparent that the results in the 
previous section are not valid. In this case, it is suggested that a Kirchhoff-Huy gen's ap- 
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Figure 4. The concentric cylindrical waveguide with 
a localized obstruction in the aperture plane x = x . 



proach is useful [Wait, 1962c]. 

It is now imagined that the incident field results from an equivalent line magnetic source 
at x=0 (i.e., 0=0). Thus, 

H(x,y)=^a m ^(t m) y)e- ixt m for x<x Q (56) 

m 

where a m is a coefficient which does not depend on x or y. It is assumed that the aperture 
plane x=x is obstructed in such a manner that the effective aperture is a slit extending from 
2/ =2/i to y 2 (i.e., r— a=z x to z 2 ). The situation is illustrated in figure 4. Thus, within the 
Kirchhoff approximation, 

H(x , y)=^a„&(t m , y)e- ix o i m for y t <y<y 2 
=0 for 0<y< Vl 
=0 for y 2 <y<y . (57) 

In other words, we are assuming that the field within the aperture of the slit has the same value 
as if the slit were not present. It is known from a study of the rigorous solutions of diffraction 
by slits that this is an excellent approximation [Born and Wolf, 1959] provided the width of 
the slit is greater than about a wavelength. 

The field in the region x^>x can now be expressed in the form 

H(x 9 y)=^2J2 A n m) Ht n , y)e- i( *-*»> t ne- 1x o<m (58) 

m n 

where 

rv2 
9A $(t m , y)$(t n , y)dy 

We see clearly that the incident mode of order m excites modes of order n where m and n are 
positive integers. 

It is convenient to write 

where, for m^n, 



and 



where 



A^ = [P^+Q^]a m (60) 

^=-^ife«r fl - wftw * (6i) 

^-^i^r^^-^' (62) 

°^=wtf) (63) 
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and 



G -^wzf) (64) 



and 



In obtaining the above forms for P£ m) and Q ( n m) use has been made of the orthogonality 
condition given by (32). 
When m=n, we have 

AP=[P™ + QP]a n (65) 

where 

P« =1 _^» [OMfdy (66) 

2/o Jo 

Qn n) ^~ 2 -^ ( V °[Gn(y)] 2 dy (67) 

where use has been made of (39). 

The integrals over the range to ij\ in the preceding equations can be regarded as the 
influence of the obstacle on the ground, whereas the integrals over y 2 to y are related to the 
protuberance at the ionosphere. To evaluate these integrals it is desirable to expand the G 
functions as a power series in ?/. 

Since 

6 n (0) = l (68) 

[dG n (y)ldy] y = Q =-q (69) 

and 

-(tn-y)G(.y)ioTmyy, (70) 



dy 
it is not difficult to show that 



Thus, 



2 3 

O n (y)G m (y)=l-2( L y+(t tl +t m +2g 2 ) |- (1+2*^+2^) f +• •• (72) 



and the expansion for [G n (y)] 2 is obtained by simply replacing t m by t n in the preceding result. 
The integrations for the P integrals are now readily carried out. They yield 



where 
and 



p ™=~j; 9m,n\yi-qy\+{t n +t m +2<f) ^-(1+2^+2^2) f§+. . .] (73) 

_ $(£ m , G) = w 2 (t n )—qw 2 (t n ) 
9m ' n *(t n ,0) w' 2 (t m )-qw 2 (t n ) {/V 

PP=l-^\yi-<Bfi+(tu+& f-(l+4^g)^+. . .]• (75) 

The Q integrals are evaluated in a very similar manner. Thus 
Q%»=~T 9 m ,nG m (y )G n (y ) [(y -y 2 ) -aidfe-lfe) 1 

+(*„+* m +2 2 ?) M--(l+2« Wi +2^ i ) ^=^+ . . .] (76) 
and 

Q r=1 _2^ [Gn(yoW [^ (?/o _ 2/2) _ gi(yo _ y2)2+(<K+g?) fa^!_ (1+4Mj) fa/?)_ 4 + . , J. (77) 
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Due to the typically large values of q t at VLF, the preceding series for the Q functions are 
probably not useful. It would be better to work directly with eqs (62) and (67). 

6. Discussion of Formulas 

Some of the previous results are now discussed briefly. For purposes of illustration, it 
is assumed that the ionosphere is a sharply bounded medium whose effective conductivity is 
€ co, where e =8.854X10 -12 . Under this condition extensive numerical values of the coeffi- 
cients t n satisfying eq (27) are available [Spies and Wait, 1961]. Using these values, the various 
quantities entering into the formulas for the modal coefficients can be evaluated in a straight- 
forward manner. 

It is seen that the modal coefficients, given by eqs (73), (75), (76), and (77) all contain 
the factor A n . This factor is a modal excitation factor and it is a measure of the efficiency of 
excitation of a given mode from a line or dipole source [Wait, 1961, 1962b]. In the present 
context it is normalized so that it approaches unity for perfect ground conductivity (g=0) 
and a flat earth (a= <» ). In general it is a complex quantity. To illustrate its behavior co r is 
set equal to 2 X 10 5 and h is taken as 70 km. Furthermore, the ground is assumed to be perfectly 
conducting. Under these conditions A n for n=l has the following complex values for the 
frequencies indicated: 

A! = 0.95 [3.8° (10 kc/s), 0.79 [3.0° (15 kc/s), 
0.59 [4.6° (20 kc/s), 0.37 [7.4° (25 kc/s), 
0.19 [10.6° (30 kc/s). (78) 

This mode corresponds to the mode of least attenuation. It is characterized by an excitation 
factor which decreases approximately as the inverse of the frequency. Under the same 
conditions A nj for n greater than 1, is roughly unity over the same frequency range [Wait, 
1962b]. 

The modal coefficient P^ defined by (75), in the case of n=l, can be written 

pa)^ 1 _2A 1 ^. (79) 

Here, hi/h is the ratio of the heights of the obstacle on the ground to the height of the ionosphere. 
This quantity would never be greater than about 0.05 and thus the modification of the first 
mode by even an extremely high mountain range would be small. This is particularly the 
case at the upper end of the VLF band where the excitation factor is small. 

The relative conversion of the field from an incident mode of order 1 to a mode of order 2 
is obtained from the factor P^ w) defined by (73) for m=l and n=2. Approximately, this 

can be written 

P^^-2A 2 g lt2 (h lt lh). (80) 

The complex quantity g lt2 is defined by eq (74) for ra = l, n=2, and #=0. For the same 
conditions, its magnitude for the frequencies indicated, are given as follows 

\g h 2 | = 1.72 (10 kc/s), 1.47 (15 kc/s), 1.08 (20 kc/s), 0.70 (25 kc/s), and 0.41 (30 kc/s). (81) 

Since |A 2 | is of the order of unity, it is thus apparent that the conversion to higher modes may be 
significant. 

The influence of the protuberance on the upper boundary is described by eqs (76) and (77) . 
The situation is similar to that of the ground obstacle except that the factors G n (y ) appear. 
Actually, these are the ratio of the field just below the ionosphere reflecting layer to the field 
at the ground. 

Of particular interest is the possibility that, as a consequence of an ionospheric irregularity, 
a mode of order 1 may be excited by an incident mode of order 2. The magnitude of this 
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first-order mode relative to amplitude of the second-order mode is obtained from eq (76) 
with m=2 and n=l. Thus, approximately 

<?f a -2A!ft. ,ft(jfe) ft (2/0) (^jp) (82) 

The numerical magnitude of the excitation factor Ai has already been discussed. As noted, 
it may be quite small for frequencies of the order of 25 kc/s. However, in certain cases, the 
height-gain function Gi(y), is somewhat greater than unity. Thus, the conversion to the 
lower-order mode may be quite significant. This point can also be demonstrated directly 
from eq (62) which, in this special case, has the form 






(83) 



Modes of the "whispering gallery' 7 type [Budden and Martin, 1962], also known as "(Mirth- 
detached modes" [Wait, 1962b], are associated with a low excitation efficiency (i.e., A t is 
small). However, the height -gain function Gi(y) for a "whispering gallery mode" is an in- 
creasing ('unction of height. Thus, the product of the integral over ?/, to //„ and the excitation 
factor Ai may be of appreciable magnitude. 

Further work on this subject awaits the completion of the extensive tabulations of nu- 
merical values of the coefficients t n and the height-gain functions G n (y) for a range of the 
parameters. Also, the complications resulting from the earth's magnetic field arc incorporated 
following the approach introduced in a previous paper [Wait, 1963]. 



I thank Prof. S. Maley and Mr. K. P. Spies for their helpful remarks, and Mrs. Eileen 
Bracket! for her assistance in preparing the manuscript. 
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